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I study the lattice dynamics and electron-phonon coupling in non-centrosymmetric quasi-one- 
dimensional K2Cr3As3 using density functional theory based first principles calculations. The 
phonon dispersions show stable phonons without any soft-mode behavior. They also exhibit features 
that point to a strong interaction of K atoms with the lattice. I find that the calculated Eliashberg 
spectral function shows a large enhancement around 50 cm“^. The phonon modes that show large 
coupling involve in-plane motions of all three species of atoms. The q dependent electron-phonon 
coupling decreases strongly away from the Qz = 0 plane. The total electron-phonon coupling is 
large with a value of Aep = 3.0, which readily explains the experimentally observed large mass 
enhancement. 

PACS numbers: 74.25.Kc 


I. INTRODUCTION 

Almost all known superconductors are conventional 
electron-phonon superconductors, and the phenomenon 
of unconventional superconductivity, especially in one¬ 
dimensional (ID) materials, is relatively unexplored. 
Therefore, the recent report of superconductivity 
that is likely of an unconventional nature in non¬ 
centrosymmetric quasi-lD K2Cr3As3 with a of 6.1 K 
by Bao et al. is of significant importance^. 

K2Cr3As3 is structurally similar to the MM03A3 or 
M2M06A6 {M = metals, X = chalcogens) family of 
compounds first discovered in 19802*^. In addition to 
K2Cr3As3, superconductivity has also been observed 
in two other isostructural compounds Rb2Cr3As3 and 
Cs2Cr3As3 with T^s of 4.8 and 2.2 K, respectively^i^— . 
The one-dimensional element in A2Cr3As3 {A = K, Rb, 
Cs) compounds is a [(Cr3As3)^“]oo chain, which is ar¬ 
ranged in a hexagonal lattice. This chain has the struc¬ 
ture of a double-walled subnanotube. The inner wall 
is made of stacked Crg octahedra that are face-sharing 
along the c axis and the outer wall is made of simi¬ 
larly stacked face-sharing Ase octahedra. The individual 
chains are separated by metal cations that provide the 
compensating charge and stabilize the structure. 

The superconductivity in K2Cr3As3 exhibits various 
features that are suggestive of an unconventional nature. 
The superconducting state arises out of a normal state 
that has a large specific heat coefficient of 70-75 mJ 
mol“^, which indicates strong quasiparticle mass 
renormalizationii^. The upper critical field exceeds the 
one-band BCS estimate of the Pauli limit by a factor 
of 3-4ii2p— . NMR experiments show that spin fluctua¬ 
tions grow as the temperature approaches Below 

Tc, the nuclear spin-lattice relaxation rate does not show 
Hebel-Slichter coherence peak that is characteristic of an 
isotropic s-wave superconductor, and the relaxation rate 
decreases with a power-law behavior suggesting the pres¬ 
ence of zero-energy excitations in the superconducting 


state. The presence of zero-energy excitations is also 
seen in the temperature dependence of the penetration 
depth, which decreases linearly at temperatures well be¬ 
low A fit of the angular dependence of the super¬ 

conducting gap to the superfluid density obtained from 
the penetration depth shows reasonable agreement with 
various p-, d- and /-wave models that exhibit nodes. 
The temperature dependence of the superfluid density 
obtained from transverse-field /rSR measurements show 
equally good fit to both an isotropic s-wave model and 
a d-wave model with line nodesi2. In addition, zero-field 
/iSR measurements show a presence of internal magnetic 
field in the superconducting state, although the magnetic 
field is more than 100 times smaller than that observed 
in Sr2Ru04. Raman scattering study of phonons with 
frequency greater than 100 cm“^ also show very weak 
electron-phonon coupling that is too small to account for 
the observed superconductivity in this materiaU2. In ad¬ 
dition to these experimental observations, various first 
principles and Hubbard model based theoretical stud¬ 
ies also find that a triplet superconducting state me¬ 
diated by magnetic fluctuations to be likely present in 
KaCraAs^ii-is. 

The experimental results that have so far been gath¬ 
ered on K2Cr3As3 rule out an isotropic s-wave super¬ 
conducting state in the weak-coupling limit, but they do 
not necessarily imply an unconventional ID supercon¬ 
ductivity mediated by magnetic fluctuations. Indeed, as 
pointed out by Balakirev et al., the superconductivity 
in K2Cr3As3 is inconsistent with such an unconventional 
superconducting stated. They find that the anisotropy of 
the upper critical field ')h{T) = H^/Hc2 is modest with 
1h{Tc) ~ 0.35 at Tc and 7 h(0) ss 1.5 at 0.7 K, which sug¬ 
gests a rather three-dimensional nature of superconduc¬ 
tivity. Even though K2Cr3As3 is structurally very one¬ 
dimensional, first principles calculations show that the 
electronic structure is decidedly three-dimensional with 
the presence of a large three-dimensional Fermi sheet in 
addition to two one-dimensional sheetsi^d^. Further- 
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more, the upper critial field parallel to the chain 
is limited by Pauli pair breaking and is indicative of 
a singlet superconductivity, whereas the upper critical 
field perpendicular to the chain shows no Pauli pair 
breaking effects and is limited by orbital pair breaking^. 
In a quasi-lD superconductivity, would have been 
limited by orbital and Pauli pair breaking effects yield¬ 
ing More importantly, superconductivity 

in K 2 Cr 3 As 3 has been observed in samples with resid¬ 
ual resistivity ratio ranging from 10 to SOii^r— , but the 
Tc is essentially insensitive to the non-magnetic impuri¬ 
ties. This suggest that the pairing is of s-wave type since 
scattering with non-magnetic impurities tends to average 
out the superconducting gap function, which would be 
detrimental to unconventional superconductors as their 
gap function has different signs in different regions of the 
Brillouin zone. 

Motivated by these observations, I have investigated 
the lattice dynamics and electron-phonon coupling in 
K 2 As 3 Cr 3 using first principles density functional per¬ 
turbation theory calculations. The phonon dispersions 
show the structure is stable, and there is an absence of 
a soft phonon behavior. Instead, they exhibit features 
that point to a strong coupling of K atoms with the lat¬ 
tice. The Eliashberg electron-phonon spectral function 
is strongly enhanced around 50 cm“^. These phonons 
correspond to vibrations of all three constituent species 
and the coupling cannot be attributed to an Einstein- 
like rattling vibration. The phonon modes that exhibit 
large coupling are characterized by in-plane motion of 
the atoms. The q dependent electron-phonon coupling is 
highly anisotropic. In particular, the coupling strongly 
decreases as one moves away from the qz = 0 plane. The 
calculated total electron-phonon coupling Aep = 3.0 is 
large and readily explains the experimentally observed 
large mass enhancement. 


II. COMPUTATIONAL APPROACH 

The phonon dispersions and electron-phonon couplings 
presented here were obtained using density functional 
perturbation theory within the generalized gradient ap¬ 
proximation (GGA) as implemented in the Quantum 
ESPRESSO packagers. I used the ultrasoft pseudopoten¬ 
tials generated by dal Gorso^i and cutoffs of 50 and 500 
Ry for basis-set and charge-density expansions, respec¬ 
tively. A 4 X 4 X 8 grid was used for the Brillouin zone 
integration in the self-consistent calculations. The dy¬ 
namical matrices were calculated on a 4 x 4 x 8 grid and 
phonon dispersions and density of states are then ob¬ 
tained from Fourier interpolation. A denser 8 x 8 x 16 
grid and a smearing of 0.008 Ry was used in the calcu¬ 
lation of the electron-phonon coupling. The electronic 
structure and structure relaxation calculations were also 
performed using the generalized full-potential method as 
implemented in the wiEN2k package^^ as a check. 

I used the experimental lattice parameters a = 9.9832 


TABLE I: Relaxed internal parameters of K 2 Cr 3 As 3 using 
non-spin-polarized pseudopotential and full-potential calcu¬ 
lations. The atomic coordinates in the P6m2 space group 
are: Asl {x,—x,0), As2 (a;, —x,0.5), Crl (a;, —0.5), Cr2 
{x,-x,0), K1 (a:,-2;,0.5), K2 (1/3, 2/3,0). 



Wyckoff 

position 

X 

pseudopotential 

X 

full-potential 

Asl 

3 / 

0.8345 

0.8344 

As2 

3k 

0.1667 

0.1668 

Grl 

3k 

0.9175 

0.9167 

Gr2 

3 / 

0.0859 

0.0864 

K1 

3k 

0.5380 

0.5381 



FIG. 1: Calculated GGA band structure of non-spin- 

polarized KaCrsAss in the relaxed structure. 

and c = 4.2304 A but relaxed the internal atomic posi¬ 
tion parameters of the P6m2 space group. The results 
of the geometry optimization calculations performed us¬ 
ing both the pseudopotential and full-potential methods 
are summarized in Table ID These are in good agreement 
with each other and to the experimental valuesi*^. Struc¬ 
tural optimization leaves the As-As distances basically 
unchanged but reduces the Cr-Gr distances by 0.1 A rela¬ 
tive to the experimental values, in accord with a previous 
theoretical studyi^. 

III. RESULTS AND DISCUSSION 

The band structure and Fermi surface of K 2 As 3 Gr 3 cal¬ 
culated using the relaxed internal parameters are shown 
in Figs. [T] and [H respectively. These differ in details 
from the ones presented in Refs. [l3 and [H that used 
experimental atomic positions. The one-dimensional na¬ 
ture of the material is evident from the relatively small 
dispersion along the in-plane T-K-M-T and A-H-L-A 
paths. However, some bands do show moderate disper¬ 
sion along these paths, which indicates that the bonding 
between [(GrsAss/^^Joo chains is substantial. The bands 
disperse strongly in the kz direction (the T-A path), and 
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FIG. 2: (Color online) Calculated GCA Fermi surface of 
non-spin-polarized K 2 Cr 3 As 3 in the relaxed structure. 


these reflect the short inter-atomic distances and strong 
bonding in the out-of-plane direction. I find substantial 
As p character above the Fermi level and the presence of 
large amount Cr d character in the lower part of the va¬ 
lence band manifold, which suggests that there is a strong 
covalency between Cr and As atoms. Even though the 
valence bands disperse relatively little in the in-plane di¬ 
rections, the valence band manifold extends from —4.8 
to 0.2 eV and has a band width of 5.0 eV. This again 
reflects the strong covalency between Cr and As atoms. 
I obtain N{Ef) = 6.94 eV“^ per formula unit both spin 
basis for a value of the electronic density of states at the 
Fermi level. This corresponds to a calculated Sommerfeld 
coefficient of 7 = 16.36 mJ/mol K^. The experimental 
7 obtained from specific heat capacity measurements is 
^70 mJ/mol K^—il, and this corresponds to an ehance- 
ment by a factor of 4.3 relative to the calculated value. 

There are three bands that cross the Fermi level, and 
these give rise to three Fermi sheets. The two upper-lying 
bands cross the Fermi level only along the kz direction 
and give rise to a pair of ID sheets. In my calculations 
these are degenerate, but they are non-degenerate in the 
calculations that use experimental atomic positiona^^d^. 
The lower-lying band crosses the Fermi level around the 
K point, and this gives rise to a large three-dimensional 
sheet. Therefore, even though the structure of K 2 Cr 3 As 3 
is quasi-lD, the material is very three-dimensional from 
the electronic structure point of view. 

The phonon dispersions and density of states (pDOS) 
of K 2 Cr 3 As 3 are shown in Figs. |3] and HI respectively. 
These again reflect the quasi-lD nature of the material 
as the in-plane dispersion of the optical branches is com¬ 
paratively smaller than the out-of-plane dispersion. The 
strong bonding along the c axis is also seen in the large 
dispersion of the out-of-plane polarized acoustic mode 
along the qz direction. This branch reaches more than 



FIG. 3: Calculated GGA phonon dispersions of non-spin- 
polarized K 2 Cr 3 As 3 . 



FIG. 4: (Color online) Calculated GGA phonon density of 
states of non-spin-polarized K 2 Cr 3 As 3 . 


100 cm“^. Although the acoustic branches have smaller 
dispersion along the in-plane direction, these do not ex¬ 
hibit soft-mode behavior. This indicates weaker bonding 
between the [(Cr 3 As 3 )^“]oo chains, but it also implies 
that the structure is stable and not in a very close vicin¬ 
ity of a structural phase transition. 

The pDOS exhibits many narrow peaks that are caused 
by narrow band widths of the phonon branches. The par¬ 
tial pDOS shows that the K atoms contribute to phonon 
modes up to 200 cm“^. Therefore, the K atoms do not 
behave like a rattler that is weakly coupled to the lattice. 
The interaction with K ions is also seen in the stronger 
dispersion around the K point for branches below 200 
cm“^. The Cr and As atoms contribute to vibrations 
all across the phonon dispersion range. In particular, 
both Cr and As atoms contribute to the high-frequency 
phonon modes that again indicates strong covalency be¬ 
tween these two. 

The strength of the interaction between electrons and 
phonons is described in terms of the Eliashberg spectral 





































































4 


function: 

'' k.q,i/,n,m 

( 1 ) 

where N{Ep) is the electronic density of states at the 
Fermi level, ej^ is the electronic energy at wavevector k 
and band index n, Wi/q is the energy of a phonon with 
wavevector q and branch index v, and 5^ k+q 
trix element for an electron in the state |nk) scattering 
to Imk + q) through a phonon w^yq. 

The calculated Eliashberg function is shown in Fig. [5] 
It shows that phonon at all frequencies contribute to 
the coupling, although the the phonon modes around 50 
cm“^ make an especially large contribution. The pDOS 
also shows a peak at this frequency. Similar low-energy 
peaks has also been observed in M 2 MoqXq compounds, 
and they have been attributed to the Einstein-like vibra¬ 
tions of M atoms^^i^i. However, we can see from Fig. 0] 
that the motion of all three types of atoms contribute 
to this low-energy peak in K 2 Cr 3 As 3 . Therefore, even if 
these phonons have relatively small dispersion, the cou¬ 
pling cannot be ascribed to an Einstein-type rattling vi¬ 
bration of K atoms. I also investigated the polarization 
of the phonons that show large coupling to the electrons 
and find that these show strong in-plane character. 

The q dependent total electron-phonon coupling Aq 
(= /p da;a^F(a;, q)/a;) is plotted for the Qz = 0 and ^ 
planes of the Brillouin zone in Fig. | 6 l I find that the 
variance of the coupling is large within a Qz plane. Fur¬ 
thermore, the coupling strongly decreases as one moves 
away from the Brillouin zone center along the qz direc¬ 
tion. This shows that the coupling is highly anisotropic, 
and it should also result in an anisotropic superconduct¬ 
ing gap function. To find out which Fermi sheets con¬ 
tribute more to the total electron-phonon coupling, I cal¬ 
culated the partial electron-phonon couplings at several 
q points by summing over the ID and 3D Fermi sheets 
separately in Eq. [T] I find that the 3D sheet makes a 
dominant contribution to the coupling. Therefore, con¬ 
trolling the size of the 3D sheet by chemical doping, pres¬ 
sure or strain might be an effective way of enhancing or 
suppressing superconductivity. It will be interesting to 
see if magnetism or charge density wave order appears 
due to nesting between ID sheets when the 3D sheet is 
suppressed. 

I obtain a large value of Aep = 3.0 for the total 
electron-phonon coupling constant Aep = X]q ^ = 

2 da;a^F(w)/a;. The logarithmically averaged 

phonon frequency is win = 71 cm“^. These can be used 
to estimate the by using the simplified Alien-Dynes 
formula, even though the use of this formula is invalid 
in the present case of strong electron-phonon coupling 
that is highly anisotropio^^. The simplified Alien-Dynes 
formula has the form 

T = exp [ ] . 

1.2 \ Aep — ^*(1 -b 0.62Aep) / 



Q) (cm'^) 

FIG. 5: Calculated Eliashberg spectral function of 

KaCrgAss. 


Qz = 0 - y 

FIG. 6 : Galculated q dependent electron-phonon coupling 
Aq shown for qz = 0 and Qz = 7 planes. The radius of the 
circle is proportional to the magnitude of Aq. 

Using a value of = 0.12 for the Coulomb pseudopo¬ 
tential parameter, I obtain Tc = 17.8 K, which strongly 
overestimates the experimental Tc of 6 K. A full solution 
of the Migdal-Eliashberg equations, which is beyond the 
scope of this paper, will be required to obtain a more 
accurate Tc based on the calculated spectral function. 
Furthermore, there might be additional reasons for the 
overestimation of Tc. In the phonon calculations, I have 
neglected the phonon anharmonicities. This can raise the 
phonon frequencies, which will lower the total electron- 
phonon coupling. Another reason for the overestimation 
might be because I have neglected the effects of mag¬ 
netic fluctuations. First principles calculations show the 
presence of various magnetic interactions with in this 
materiali^ii^, and this will also suppress Tc. 

The calculated total electron-phonon coupling Aep = 
3.0 readily explains the experimentally observed large 
electron mass enhancement from specific heat measur- 
ments. The electron-phonon coupling enhances the mass 
by a factor of 1 -b Aep = 4.0, which is very close to 
the experimentally observed mass enhancement factor 
of 4.3. Therefore, the large mass enhancement in this 
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material does not come from strong electronic correla¬ 
tions (in the sense of electron-electron interactions) or 
spin fluctuations. This is consistent with the experi¬ 
mentally estimated Wilson ration of Rw ~ 1— that sug¬ 
gests a weak magnetic enhancement. The strong cou¬ 
pling nature of superconductivity in K 2 Cr 3 As 3 is also 
evident in the large jump of the dimensionless specific 
heat AC/yTc = 2.4^*^ at Tc- The strong electron-phonon 
coupling superconductivity proposed here can further be 
confirmed from tunneling experiments that should show 
a large value for the ratio 2A(0)/Tc, where A(0) is the 
superconducting gap. 

Although the calculations presented here imply that 
the superconductivity in K 2 Cr 3 As 3 is caused by electron- 
phonon interactions, this does not necessarily mean that 
the superconductivity is dull or of the conventional BCS 
type. The crystal structure of this material lacks inver¬ 
sion symmetry, and parity is not a good quantum number 
is such a case. The pairing will likely have a dominant s 
character mixed with some p character. This can explain 
the presence of weak internal magnetic field below in 
/iSR measurements^^. In addition, the electron-phonon 
coupling is highly anisotropic, which can lead to acciden¬ 
tal gaps in the superconducting state. This would be con¬ 
sistent with the signatures of nodes in the gap function 
in the NMRpi^ and penetration depthii measurements as 
well as the absence of Hebel-Slichter peak in the NMR 
measurements. 

The presence of such a large electron-phonon coupling 
in this material provides an opportunity to study the 
properties of a superconducting state in the strong cou¬ 
pling limit. In particular, this may help explain why 
the measured upper critical fields and h }.2 greatly 
exceeds the Pauli limit estimated for a weakly-coupled 
isotropic BCS superconductor. It would be interesting 
to see if the fact that the coupling strongly decreases 
away from the (/z = 0 plane can explain why 
in this material. These issues require more theoretical 
studies based on the full solution of Migdal-Eliashberg 
equations to be fully clarified. 

As mentioned above, superconductivity is also ob¬ 
served in other inorganic quasi-ID compounds includ¬ 
ing Tl 2 Mo 6 Se 6 and Lio.gMoeOiy, and the nature of their 
pairing interaction has not been fully elucidated. The 
results of phonon dispersions and electron-phonon cou¬ 
pling for K 2 Cr 3 As 3 presented here provides a reference 


case for other quasi-lD superconductors. In particular, 
the present study shows that a behavior that deviates 
from what is expected from a weakly-coupled isotropic 
s-wave BCS superconductor does not necessarily imply 
that the superconductivity is mediated by magnetic in¬ 
teractions. 


IV. SUMMARY AND CONCLUSIONS 

In summary, I have investigated the lattice dynam¬ 
ics and electron-phonon coupling in K 2 Cr 3 As 3 using first 
principles calculations. The phonon dispersions do not 
show closeness to a lattice instability and exhibits fea¬ 
tures pointing to a strong coupling of K atoms with the 
lattice. The Eliashberg spectral function is strongly en¬ 
hance around 50 cm“^, but there is coupling to high fre¬ 
quency phonons as well. The peak in the spectral func¬ 
tion corresponds to the vibrations of all three constituent 
atomic species of the material. Although the electrons 
do not predominantly couple to vibrations of any specific 
atoms, vibrations that involve in-plane motion of atoms 
show strong electron-phonon coupling. The electron- 
phonon coupling shows strong q dependence, especially 
along the out-of-plane direction. The calculated total- 
electron phonon coupling Aep = 3.0 is large and is con¬ 
sistent with the experimentally observed large mass en¬ 
hancement and jump in the specific heat at Tj^. These 
calculations motivate tunneling experiments that can fur¬ 
ther confirm the strong coupling nature of superconduc¬ 
tivity in this material by displaying an enhanced value of 
the ratio 2A(0)/rc- The electron-phonon spectral func¬ 
tion can also be extracted from such experiments, and 
they can be compared with the results presented here 
to validate the strong electron-phonon coupling scenario 
proposed here. 
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